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Abstract 



The oscillator representation method is presented and used to calculate the 
energy spectra of a superposition of Coulomb and power-law potentials and for 
Coulomb and Yukawa potentials. The method provides an efficient way to obtain 
analytic results for arbitrary values of the parameters specifying the above-type po- 
tentials. The calculated energies of the ground and excited states of the quantum 
systems in question are found to comply well with the exact results. 



PACS numbers: 03.65.-w, 32.30.-r, 32.60.+i 



1 Introduction 

The combination of Coulomb potential with different type of confining potentials is a typ- 
ical problem in quantum mechanics for many applications. The two-dimensional electron 
gas with Coulomb interaction confined in a parabolic potential is a model of a quantum 
dot (see for example M, 12] and references therein), i.e. an artificial atom the properties 
of which can be controlled by a man (for review on quantum dots see |3j, ||J]). The com- 
bination of Coulomb and Yukawa potentials, called the exponentially screened Coulomb 
potential (ESCP), describe the effective two-particle interaction for charged particles in 
an ionized gas or in a metal (see || and references therein). The model of quark confine- 
ment of charmed quarks is successfully simulated by a simple superposition of the linear 
and Coulomb term M. 

Whereas many approximate analytical and numerical methods have been developed 
for the last years for the analysis of such a superposition of potentials [0, we would like 
to suggest a method which is called the oscillator representation method (ORM) |§. The 
ORM, based on the ideas and methods of quantum field theory, has been proposed to 
calculate the binding energy of different systems with fairly arbitrary potentials described 
by the Schrodinger equation (SE). A spectrum of the one-dimensional power-law potential 
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or of the quartic oscillator, which are test benches for approximate methods (see for 
example || |TD[|) , can be determined analytically in the ORM with accuracy very closed 
to the best, available nowadays, computational methods (see for the comparision j8|). 
One of our goals is to present the mathematical formalism of the ORM in detail making 
it available for further applications, in particular, for calculations of ground and excited 
states of the systems with Coulomb interaction. 

In situation, when different terms of the effective potential of a system can compete, 
and there is no basis for the perturbative approach to be applied, a solution of the SE 
becomes a nontrivial problem. However, the lowest levels of such a system, especially ones 
closest to the ground state, can be described with a good approximation by the oscillator 
states defined in the expanded hyperspace which includes the original one. For the first 
time, Fock demonstrated this result for a nonrelativistic hydrogen atom by making a 



stereographic projection of the momentum space onto a four-dimensional unit sphere [11]. 
For an arbitrary smooth potential, admitting the existence of system bound states, its 
eigenfunctions have exponential asymptotic behavior at large distances JT2], [13| . Therefore, 
the variables in the original SE must be changed so that the modified equation should 
have solutions with the oscillator behavior at large distances. Schrodinger pointed out 
at the existence of such a transformation which turned the three-dimensional Coulomb 
system into the oscillator one in the four-dimensional space |14J]. The explicit form of 



this transformation has been found in |T5[ and used to solve the classical Kepler problem. 
Since this transformation is not a canonical one, after the transformation we have a new 
system with another set of quantum numbers and wave functions which contains, however, 
a subset of the original wave functions. The transformation of variables, leading to the 
Gaussian asymptotic behavior in the expanded space, is one of the basic elements of the 
ORM. 

As the next step, it is necessary to represent the canonical variables (coordinate and 
momentum) of the Hamiltonian through the creation and annihilation operators a + and 
a. The pure oscillator part with some, yet unknown, frequency u is extracted from the 
Hamiltonian, i.e. H — > Hq + Hj = ua + a + higher order terms. The remaining part, 
i.e. the interaction Hamiltonian Hi, is represented in terms of normal products over a + 
and a. In addition, it is required that the interaction Hamiltonian does not contain terms 
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quadratic in the canonical variables. This condition is equivalent to the equation 

^ = (1) 
duo 

which determines u, the oscillator frequency, in the ORM and is called the oscillator rep- 
resentation condition (ORC) ||. Similar ideas are used in the Hartee-Fock-Bogoliubov 
theory to describe different correlations between nucleons moving in an average nuclear 
potential (see for review |I5|), in order to include particle correlations into a new quasi- 



particle vacuum. In fact, the ORC coincides with the variational equation arising after 
averaging of the total Hamiltonian over the Gaussian probe functions in the expanded con- 
figuration space. In addition, the dimension of the hyperspace is the variational parameter 
in the ORM. Using the ORM, we will calculate the energy spectrum and wave functions for 
Coulomb and power-law potentials and for the exponentially screened Coulomb potential. 

The paper is organized as follows: In Sect. 2, basic formulae of the ORM for spherical 
potentials are presented. In Sect. 3 the method is applied to calculate the energy spectrum 
for Coulomb and power-law potentials, for Coulomb and Yukawa potentials. In Sect. 4 
we discuss the main results. In the appendices, useful relations used in the ORM are 
presented. 

2 ORM in the Space R D 

Let us consider the Schrodinger equation in 3-dimensions 

¥(r) = (2) 



J d 3 r^{r) 



— -A + V(r) — E 
2 



where V(r) is a potential of Coulomb and/or Yukawa type, decreasing at large distances 
(V(r) ^ as r — > oo). If ^(r) = ip n i(r)Y n i(9,<p), this equation for the wave function 
of the l-th orbital excitation is |12|, |13| 



f°° , / , / xn r 1 / d \2 1(1 + 1) XT , s 
I dr{r^ nl {r)) --(-) + -k_i + V( r ) - E 



{ri>ra(r)) = (3) 



where the function if) n i(r) depends on a radial variable only. General discussion on the 
transition to arbitrary dimensions in the radial SE is given in |17]]. In this case, calculation 



of the wave function ipoi(r) would be equivalent to that of the ground state wave function 
of the modified Hamiltonian in another dimension. The radial excitation wave functions 
ipni(r) — \n) (n = n r + I + 1) are equivalent to the highest oscillator states. 
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At large distances, for Coulomb and Yukawa type potentials the asymptotic behaviour 
of the wave function ij){r) as r —>■ oo is well known, i.e. ip(r) ~ exp(— a(r)) . The new 
variable r = r(q) should be introduced so that a(r(q)) ~ q 2 as q — > oo and 
ip{r(q)) ~ exp(— q 2 ). Using the substitutions 



r = q 



2p 



and rip i(r) = q a (p(q), 



(4) 



after simple calculations we obtain 



dqq D - l 4>{q) 



1/ / d 



2 V \dq 



D- 1 d 
q dq 



Wj(q 2 -E) 0(g) = 



where 



(5) 

(6) 
(7) 



W^E) = +Ap 2 (q 2 Y 2 ^\V(q 2 ) - E), 

D = 2a - 2p + 2 

K(l,p,D) = ±((D-2) 2 -4p 2 (2l + l) 2 

We can identify Eq.([|) with the SE in the expanded space R D for the wave function <j)(q) 
depending on the radial variable q only 



rf D q0(q) 



-A D + W I (q 2 ] E)-e(E) 



(q) = o, 



where the function 



e(E)=e(l,p,D;E) 
should be considered as an eigenvalue of the SE in D-dimensions 



-A D + W I (q 2 ,E)}<j ) (q)=e(E) ( f ) (q) 



(9) 



(10) 



We note that the potential Wj(q 2 ; E) contains the attractive term which can compensate 
the repulsive part of the V(q 2 ). 

The energy spectrum E n i of the original system is contained in the radial excitation 
spectrum e' 71 ^ of the Hamiltonian of Eq.flTTD 

E{E)^\q) = e^\E)^\q) , (n r = 0,1,2,...). (12) 



and it is determined by the equation 



e [nr \E) =e [nr] (i,p,D;E) = 0. 



(13) 



The additional parameters p and D can be found, for example, by the minimization 
of the energy in the zeroth approximation 

Eq(E) = min e(l, p, D; E). 
{p,d} 

On the other hand, the parameter p is connected with the condition of the Gaussian 
asymptotic behavior for the wave function (see the discussion above). The parameter D 
is connected with the behaviour of the wave function at short distances. If the potential 
V(r) has no a repulsive character as r — > in the ORM we choose K(l, p, D) = and 

D = 2 + 2p(2/+ 1). (14) 

In the opposite case, the parameter D is chosen to suppress a repulsion produced by 
the potential V(r) at r — > 0. In fact, the parameter D is used to improve the zeroth 
approximation at the minimization of the ground state energy. 

One can see that the radial quantum number n r does not enter into Eq. (|TT|) in an 
explicit form. The orbital quantum number i is absorbed by the parameter D. From the 
point of view of the expanded space R D the functions 

<t>n(q) = q 2p - a ^m(q 2p ) or Vw(r) = r 2 ^^^) 

for any n and for fixed I are eigenfunctions of the basic series of radial excitations in the 
space R, D with radial quantum number n and zeroth orbital momentum. Consequently, 
the solution of the SE in 3 dimensions for I orbital excitation is equivalent to the solution 
of the SE in the space R D for states with the zeroth angular momentum. As a result, the 
original SE is represented in the form of Eq. (|TTD in which the wave function of the ground 
state 4>{q) has 

• a Gaussian asymptotic behaviour at large distances 0(g) ~ exp(— q 2 ), 

• a maximum at the point q = 0. 

In addition, the integration over the expanded space R D and the Laplacian are defined 
by: 

. 2tt d / 2 °r 



5 



, . ( d 2 D-l d\ 



The next step is to introduce the oscillator representation which can be done in the 
following way. The Hamiltonian H can be rewritten in the form 

H = \{p 2 + co 2 q 2 ) + (W(q) - \u?q 2 ) , (15) 

where uo is an unknown oscillator frequency. Let us substitute the creation and annihila- 
tion operators (Appendix A) into Eq. (|15|) and go to the normal product of the operators 



alj and a,j. Hereafter we use the notation d = D. One can obtain 

hj? + u 2 q 2 ) = uJ2 a j a j + ^ = uj(a + a) + ^u, (16) 

3 

where : * : is the symbol of the normal ordering and (qk) = J2j kjQji 

W d (k 2 )= f(dp) d W(p)e^ . (17) 



We require that the interaction part of the Hamiltonian should not contain the term with 
: q 2 : because this term is postulated to be included into the oscillator part completely. 
This requirement gives the equation for the frequency u: 

^_/(*)^ )eX p ( -^ = 0. (18) 

Using these equations, we can rewrite the Hamiltonian, Eq.(|l|), in the form: 

H = H + Ej + e , (19) 

H = uj(a + a) , (20) 



Hr = JQ d W d (k 2 )eM-^)--e^-l + 



2d 



V"V(-^) : exp(-g 2 + 2{pq)) - 1 + q\l - %) : . (21) 



duj f /dk\d~ - 9 . , k 2 . . . 

£ « = T + iU ^(fc 2 )exp(--). (22) 



Details of calculations of different integrals appearing in the above expressions can be 
found in Appendices. In the ORM the solution of the SE Eq.flTTD has the form 

<t> n (q) = exp(-^g 2 ) ^ c^P^ (? V> > 

where {P^ {t)} is the class of orthogonal polynomials which are orthogonal within the 
interval < t < oo with the weight function 



X 

jdt t%- 1 e- t Pi D \t)P! n D \t) = 5 nm . 



o 



These orthogonal polynomials can be constructed by using the formalism of creation 
and annihilation operators dj and Oj in the space K D (see Appendices A, B). Using the 
definition for the ground wave function in the R D space 



10) = II = (-) ^ ( 23 ) 

we can write for radial excitations n = 1, 2, ... 



3> n ~ (a + a+) n \0) ~ Pi D \uq 2 )e-^ 2 (24) 
~ Pi D \ur 1 /p)exp(-Y 1 /p), 

where P^ D) (t) is a polynomial of the n th order. The parameter D in this representation 
can be considered to be any positive number. These polynomials satisfy the orthogonal 
condition 

($ n ,$ m ) ~ (0\(aana + a + ) m \0) 

dq g D - 1 exp(-g 2 )Pf)(g 2 )p(f)(g 2 ) 



o 

oo 



~fdt t D ^eM-t)PW(t)Pi D \t)~5, 
o 

Thus, the ORC, Eq.(P, is written as 





, ) j (E;uj,d) = . (25) 



This equation determines the parameter u = ui(E, d) as a function of the energy E, d 
and other parameters defining the potential V(r) in Eq.(|3|). 

According to Eq.(|Kj), the ground state energy E of the original problem in the iV-th 
perturbation order of the ORM is defined by Eq.(E3) and by the equation 



e (N) {E, d) = e {E, d) + e 2 (E, d) + ... + e N (E, d) = . (26) 

Eq.(|26"D determines the energy E^)(d) in the N-th perturbation order as a function of 
d or p and other parameters defining the potential. The function e(E;u,d) in Eq.fll9|) 
depends on two parameters d and u. Let us denote the parameter d and other additional 
parameters by {a,}. The ground state energy in the expanded space up to second order 
is 

£(2)(E;uj,aj) = e (E;u,aj) + e 2 (E;u,aj) . 

According to Eqs.(|22]), fll7l) and (||), the function £q{E; u>, acj) can be written in the form 

so(E; u, otj) = A(w, atj) — E ■ B(u, ctj), (27) 

where A(u,atj) and B(u,otj) are known functions. 

In the zeroth order Eqs.(P^D and ( $T\) lead to equations 

^-E (u;, a 3 ) = A (^r^\) = for all fc, (28) 
oa k da k \B{uu,oij)J 

which determine the parameters {ctj} as functions of Eq 



a 



j 



aj (E). (29) 



As a result, in the zeroth approxiamtion we have 

Ea = miD = A f>^; po) 

{u>, aj } B(ui, atj) B(uj , a)) 

where the parameters ujq and a° determine the minimum. 

In the second approximation we solve the following equation 

S{2){E] Oil, aj) = e (E; a>, acj) + e 2 (E; u, atj) = . (31) 

We assume that the second correction is small one and the energy is Em = Eq + E 2 . 
Consequently, we have 

e (E {2 y, oo , or") = A(u , a?) - E (2) B(uj , a]), (32) 
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where uq = u>(E ). Therefore, the second correction is 

E 2 

Finally, we have 



e 2 (E Q ;u ,a° j ) 
B(u ,a°) + ° {E ^ 



E(2) = Eq + E2 = 

A(aij) e 2 (E ;ujQ,a j (E )) 



(33) 



= mm 



{J,a-} B(oij) B(u ,aj(E )) 



+ E O 





E2 




Eq 



B(uj , a C j 



^ + E O 





E 2 


( 


Eq 



In our calculations we use essentially these equations. The accuracy of the oscillator 
representation can be evaluated as 

A £2 
~ — 

^0 

For the radial excitations \n r ) (n r = 1, 2, ...) the matrix element 

(n r \Hj\n r ) = A lnr] (w,aj) - EB [nr] (u,aj) ^ . 

The energy in the lowest approximation looks like 

e\ w \e) = (n r \H\n r ) = e (E) + 2n r u + (n r |i?/|n r 
= A^ r \uj, otj) - EBi^iu, a 3 ) , 



(34) 



where 



Two equations 



A^icu, aj) = A(u, olj) + 2n r uj + A [nr] (uj, aj) , 
B[ nr \uj, aj ) = B(u, aj) + BH(u, a 3 ) . 



d d 
— A(u, aj) - E—B(uj, aj) = , 



(35) 



A ( ( lr \u;,aj)-EB{ nr) (u;,aj) = , 

determine the functions uj(aj) and E(aj). The energy of the n r -th excited state in the 
first approximation of the oscillator representation is determined as 

Ki B\ lr \uo{aj),a 3 ) 
Further steps should be done according to the rules formulated above. 
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3 Coulomb type potentials 

3.1 Coulomb and power— law confining potentials 

In the general case, the potential under consideration can be written as follows 



V{r) 



— + g-r 

r 



(36) 



At v — 1 this potential corresponds to the well-known "Cornell" potential ||, used to 
study the energy spectrum of a many-quark system. In addition, at v = 1 and 2 the 
potential can be used to study the Stark problem and the Zeeman effect in a spherical 
system, respectively. Due to various applications of this kind of potential, there are many 
attempts to solve the eigenvalue problem with a good quantitative accuracy (see, for 
example 0, [18|, and references therein). In fact, the main approach is based on ideas 
of the perturbation theory. As a result, the problem of summing the perturbation series 
is not solved completely. 

We apply the ORM to calculate the energy spectrum for an arbitrary value of the 
parameter v. Using the transformation Eq.(|4]), after some simple calculations the SE 
with the potential Eq . (|36|) can be written in the following form 

1 / d 2 d - 1 d ' 

2 [dq^ 4 



4p 2 Eq 



2(2p-l) 



q dq / 

+Ap 2 gq 2{vp+2p - 1) ] $(g) = e(E)$(q) . 



4p z q 



2„2(p-l) 



(37) 



d = 2 + 2p + Apl 

According to the ORM, the Hamiltonian has the form of Eq . ([19|) where 
H = oo (a + aj , 



(38) 



(39) 



dx 



-75pl e 



Ap 2 E x 



-2p 



4p 2 x~ p Agp 2 



x 



-2p—vp 



u 2 P -i r (! _ 2 P ) ujp- 1 r(i - p) u^p- 1 r(i -2 P - up) 

and the function Eq is 



e (E) = mine (^;a;,p) 
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Here 



e (E;u,p) 



duo AEp 2 T(d/2 + 2p-l) 
T 



2p-l 



r(d/2) 



(40) 



4p 2 r(d/2 + p-l) 4#p 2 r(d/2 + 2p + z/p - 1) 



ujp- 1 T(d/2) uPp+vp- 1 r(d/2) 

The energy of the system and the new oscillator frequency are determined from Eq.flTS]) 



and from the ORC, Eq. (p5|) , respectively . After some simplifications, we obtain the 
following equation from the system of two equations: 

' Z 2 r(d/2 + l) ZT(d/2 + p-l) 



E 



mm 



+ 



[ P ] [8p 2 T(d/2 + 2p-l) T(d/2 + 2p-l) 
g T(d/2 + 2p + vp-l)\ 



Z" T(d/2 + 2/0-1) 
The parameter Z is determined from Eq. (|25|) 



(41) 



z 2+u 



Z 1+l/ 4p 2 



T{d/2 + p-l) 
T(d/2 + 2p- 1) 



(42) 



guAp 



2 T(d/2 + 2p + vp-l) 
F(d/2 + 2p - 1) 



and the oscillator frequency is 



Z l ' p . 



In the strong coupling limit, i.e. as g — > oo, we obtain from Eq.(£|l]) the following 
result 



where 



C = min 

[p] 



E = gi+» ■ C 



1 1 

2 + ~, 



T{d/2 + 1) 



4p 2 T(d/2 + 2p-l) 



4z/p : 



T(d/2 + 2p + up- 1) 



T(d/2 



2 

2+i> 



(43) 



(44) 



In Table 1 we present, as an illustration, the results at the zeroth approximation of 
the ORM for v = 1. 
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3.2 The exponentially screened Coulomb potential 

The ESCP describes various physical phenomena. It was used to evaluate the cohesive 



energy of alkali metals [ 20fl , to calculate of energy bands of sodium by empirical pseu- 



dopotential method pi] and to describe of the exciton and donor spectra in ionic crystals 
and polar semiconductors [23| . 

The Schrodinger equation for the ESCP has the following form: 



A 



2/i 



P 



+ B 



\l/(r) = E^(r) , 



(45) 



where p is relative mass, A, B and c are some constants. We use the following units: the 
length unit is h 2 /(Ap) and the energy unit is pA 2 /(2Ti 2 ). The SE in the dimensionless 
units can be written 



-Ia 

2 



1 Be 

r 



2 r 



E 
~2 



V(r) = 



Using the transformation Eq.@, we obtain the modified SE 

d-1 d 



- 2p z Eq 



2(2p-l) 



1 (cP_ 

2 \dq 2 q dq 

+2Bp 2 q 2{2p - 1] exp{-cg 2p }] = e(E)$(q) ; 
e(E) = , d = 2 + 2p + Ap£ . 

Applying the ORM to the Hamiltonian Eq. (f47|) , we obtain 



2„2(p-l) 



H = H + Hi + Eq, 



where 



Hi 



10 
oo 



dx 



-2p 2 E x 



-2p 



4p 2 x~ 



u 2 P -i r (! _ 2 p) up- 1 T(l - p) 



+2Bp 2 £ 



dx 



j Q kiuj^+p- 1 J r(i - p(i + k)) 



(46) 



(47) 



(48) 
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duo _ 2Ep 2 T(d/2 + 2p - 1) _ 4p 2 T(d/2 + p - 1) 



+ fp)l fare 1J «pt- c UJ>- (49) 



o 



The energy E and the new oscillator frequency of the system are determined from Eq. ([T3|) 
and Eq.(|25])) respectively. After some algebraic calculations, the solution of the system of 
these two equations has the form 

. f Z 2 T{2 + p + 2pl) £(2p + 2p0 

E = Z n W ns„+2pt) - 2Z r( 3p+ 2 P () (50) 

+ r(3^/*" wl ^(-"-r'> 

Here the parameter Z is determined from Eq. 



*-Vf*&±™ (51) 
r(2 + p + 2pt) 



T(2 + p + 2p£) 
and the oscillator frequency is 



z l / p . 



In Tables 2 and 3 we present the result of calculations of the ESCP energy for different 
sets of parameters. The comparision with the computational results [[J is very good and 
it gives a strong support to the ideas of the ORM. 

4 Summary 

We presented the main ideas and the mathematical formalism of the ORM which is 
very efficient in getting the analytical solution to bound states of a stationary quantum- 
mechanical problem. 

Using the ORM, the discrete eigenvalues for the superposition of Coulomb and power- 
law potentials have been calculated. The results are obtained for arbitrary values of the 
strength and the degree of the power-law potential. It is a well-known fact that many 
methods fail to solve the problem of transition from a weak to a strong coupling limit. 
The advantage of the ORM is in the consistent treatment of the problem, while different 
approaches run into a problem of summing of the perturbation series. In addition, the 
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ground state and the orbital and radial excitations are treated in a similar way within the 
ORM. We note that our results are in the remarkable agreement with the ones, obtained 
for example in [23] at different limits. The ORM has been applied successfully to calculate 



the energy spectrum of Coulomb and Yukawa potentials for different sets of their strengths. 
The results show good agreement with the quantitative solutions of the problem In 
fact, in both cases our solutions have been obtained avoiding the question of the Dyson 



phenomenon [g5j which causes the main difficulties in the perturbation approach. We 
hope that the results obtained for the considered problems can be useful for an analysis 
of different questions arising, in particular, in atomic physics. 

In conclusion, the ORM allows in many cases to obtain analytical results, which are 
available only numerically or partially analytically within the perturbation approach or 
the WKB approximation. It is clear, that any progress in the development of analytical 
methods can provide a deeper insight into solution and understanding of the quantum 
mechanical problem, making further contribution to quantum mechanics. 

A Calculation of some products of the annihilation 
and creation operators 

The canonical variables qj and pj are coordinates and momentum, respectively, in the 
d-space R d and are expressed through a+ and aj 

1 / j_\ [UJ cij — a~j 



q> = (aj + aj) , Pi = V 2" ■ (A ' 1} 

The operators a~j~ and aj fulfill the standard commutation relation 

5 itj , j = 1,2,.. A. (A.2) 



Hereafter d = D. From ( |A. 1| ) and ( |A.2| ) we can obtain for q 2 , p 2 



2^,2 2 dui 2 / a o \ 

q = — + : q : , p = — + : p : , (A.3) 

where the symbol : * : means the normal ordering and u is the oscillator frequency. 
For Oj and aj we have 

e ika e ipa+ _ e %pa+ ^ka _ g-(fep) £\ 

where k and p are vectors in the d-space. 
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Let consider the expression 

Yj(k) = A+e" ife . (A.5) 

At k = from Eg. ( |A.5| ) we have 

Yj(0) = 4- ( A - 6 ) 
Taking into account ( A. 2 ), we obtain from (|A.5|) 

d im _ i5 , (A7) 

Integrating over k\ and taking into account (|A.6|) , we have 



Yj(k) = e lM a+e- ikd = a+ + ikj . (A.8) 

Similar expressions can be obtained for 

e-^ a \e^ = aj + i Pj , (A.9) 
e aala aj e- aalci = a je - a 
e aa ' +d a+e- aa ' +d = a+e a . 

According to the ORM, we transform the Hamiltonian rewriting all variables into the 
d-space in a normal order. As a result, the potential can be expressed as a polynomial of 
different powers of q. Some typical cases of potentials are the following: 

a) a power-law potential 

<F = {-IT— ■ e^ 2 lx=0 (A.10) 



d n f ( dr\ 



a 



1 r(f + n) 2 _ n r(f+n) 



" cu" r(f) 1 ' q lu™- 1 r(f + i) 

(-If rf n r ( dr]\ d 2 {1+x) _ -2iy/w{ m ) . 



7T 



where n — 1, 2, ... are integer and positive numbers 
b) an inverse power-law potential 



o 



15 



dx 



T(-t) 



-x 



l-T 



/7F / 



1 r(f + r) 



^ r(f) 



2. r T(^ + r) 



T-l 



r(f + i) 



1 ? & 1 _ 

+ — / rX 1 T 

v y 



Here r ^ n; we use the notation ef = e x — 1 — x — x 2 /2. 
c) In the general case we use the Fourier transform 



W(q 2 ) 



'dk\d~ 



2vr 

^) d W d (k 2 )ex P Uk 
All' v 



a + cT 



/ (^)'W 2 ) ex P ( - 4j) ex P ( ik ^) exp 



— ) V d (fc 2 ) exp ( - — ) : exp(ikq) 



where (kq) = J] kjqj and 



(A.12) 



W d {k 2 



(dpfW(p 2 )e 



The normal order form of an arbitrary potential can be defined by Eq s . (|A . 1 0|) - ([A . 1 2|) . 
In particular, we have 



for n = 1, 2, 3 



d 



2 "> 1 2 
4 d(d + 2) d + 2 2 4 



(A.13) 



4cu 2 



6 d(d + 2)(d + 4) 3(d + 2)(d + 4) 2 3(d + 4) 4 , . e 



2u 



For a calculation of the energy and different matrix elements we need to know the 
result of action of a combination (a + a + ) n or (aa) n . The following representation is useful 



(a + a + ) n = (-1)' 



-I3(a+a+) 



■l) n 



dP 



d(3 r - 



d n 



(A.14) 



7T 



13=0 
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According to Eqs. ( |A. 10|) - ([A . 1 2| ) , the interaction Hamiltonian contains the term pro- 
portional to exp 2 {— i(ka)}. Let us introduce the following operators: 



-i(ka) _ p -iv(ka) 



where P v is the operator defined according to the following rules: 

P v ■ Const = 

P u ■ v n = , at n < 2 

P v ■ v n = 1 , at n > 2 



(A.15) 



(A.16) 



B Normalization 

In the ORM, the wave function is defined 

\n r ) = C nr (a+a+) r |0) , j = h d, 



(B.l) 



where C Hr is the normalization constant. Hereafter we use the notation n = n r . The 
constant can be determined from the condition 



1 =< n \n) = C£<0| (a iai ) n (a+ at)" |0) 



(B.2) 



Using the definition ( |B.1|) , the representation ( A.10 ), and the normalization condition 
< 0|0) = 1 from (|B]1) we obtain 



l = Cl 



2/i. 



dr) 



d£ 



-v 2 -e 



n da n d(3 n 
x(0|e" 2 ^ (<) ■ e - 2i ^ a+v) \0) 



71 



a,/3=0 



(B-3) 



C 



0' 



2d 



c 



n da n d(5 n 

Q2n 



dr] " 
/? 
1 



n da n df3 n (l-Aapy/ 2 
Finally, for C n we can write 

' r(d/2) 



7T 



a,/3=0 



a,/3=0 



1/2 



A n n\T(d/2 + n) 



(B.4) 
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Let consider the result of the action of some operators on the wave function (P-l| ). By 
definition ((0|a+ = 0, a|0) = 0, where (0|0) = 1), it follows that 



(B.5) 



CJ-1) 



d(3 n 



drj 



71 

/3e~ 2B 



10} 



- V 2 e -B(a+a) e -2i^/f3(a+r,) I q\ 
-2nB 



(3=0 



(3=0 



\n) 



C Calculation of corrections to the energy spectrum 



In the ORM the Hamiltonian has the form 



H = Hq + Hj + sq . 



(C.l) 



The contribution of the interaction Hamiltonian Hj is considered in the perturbation 
approach. The energy and wave function are defined 



E n = E^ + E^ + E^ + E^ + ... 
*n = ^ 0) + *i 1) + *i 2) + --- , 



(C.2) 



where = £o is the ground state energy. Within the perturbation approach the SE has 
the form 



k\+k 2 =k 



(*a) 
n i 



(C.3) 



where n is a radial quantum number and k is an order of the perturbation approach. In 
the zeroth approximation, we have from Eq.( C3 ) 



H y(0) = £(0)^(0) 
D n n n 



(C.4) 



According to the ORM, H = u(a + a) and \E^ 0) = \n >. From Eq.(|Cl|) it follows that 

E W = (*ftf *(°)) = ( n \ u (a+a)\n) , (C.5) 
H = u{a + a) , = |n) . 

For k=l we obtain from Eq.( |C.3| ) 



U 1 n 1 I ^ n n n 1 n n 



(C.6) 
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Taking into account Eq.( |C.5| ), we obtain from Eq. ( |C.6|) 

n \ n I n 



(C.7) 



H n — Inuj 



Hj - £«1 ^ (0) . 

1 n n 



Following the same procedure, from Eq. (|C.3|) we obtain 

E^ = - 



Hi - Eg*] 


1 


Hj - E^_ 






Hq — 2nu 







(C.{ 



Calculation of 

Taking into account Eq. (|B.l|) and Eq.( A.14 ), we obtain from Eq.( 



2n 



E n 0) = ^(Olia^nafa^atatriO) = wC*^^ 



(C.9) 



dr]\ [ f d£, 



71 



e -^-« 2 (0|e- 2i ^ a ^ (a+ aj )( 



a,/3=0 



Using Eqs.Q, Q) and flA~9l) , after some manipulations we have 



= coC 2 



d 2n d f ( dr] y 
■da n dp n "drJ \^/tt / 

Q2n q 1 



a,P=0; t=1 



' n da n df3 n dr (i _ 4r 2 «/?) d/2 
Finally, we obtain from Eq. flC.10l) 

= (n\H \n) = 2nto . 



a,/3=0; r=l 

(CIO) 



(c.n; 



In order to demonstrate the principles of a calculation of the first order, let us 
consider the interaction Hamiltonian in the form 



H! = g 



'-IY d T 



1 ' . e 2 



uj t dx T J \\/tt / 
where g is the interaction constant and r is some constant. 

Calculation of E^\ 

Taking into account Eq. flC.l2D , we obtain from ( |C.7| ) 



2=0 



(C.12) 



E W 



■IV d T 



dr] 



-7) 2 (l+x) 



(n\ : & 



-2i\/xu>(qn) 



: \n) 



x=0 



(C.13) 
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For calculations of the first order contributions, we very often treat the following 
term 

T n (x) = [ ( ^L) d e-^ 1+x \n\ : e 2 2V ^ (w) : \n). (C.14) 



7T 



Let us consider in more detail the procedure of calculating this term. Taking into 
account Egs- flOp , Q , ([Alii) , we obtain for Eq.( |ClD 



dr) 



-T) 2 (l+X), 



IT 



e ^->(n\ : e"^ 1 ™ : \n) 



(C.15) 



9 



,2/i 



7T 



7T 



Finally, we have 



In n 



dv\ f fd£x\ f ( d&V v 2 (1+x) _ e 2 +i 2 



a,(3=0 



l da n dp> 

^ e -2i^/a(a^ 1 ) e -iuV^(a+ri) e -iuy^(ari) e -2i^(a+^2) |Q^ 



Ux) = EEC" 1 )*?!" 



rn+n) 



, , x + x) fc + d / 2 r(n + d/2) r(n - s + 1) 
T(A; + n-s + d/2) 

r 2 (k - s + i)r(2s - fc + 1) ' 

According to Eq.( |C.13| ), the first order correction to the energy of a state is 

'-l) r d 1 



(C.16) 



T r, T -Tn{x) 

UJ T ox 



x=0 



(C.17) 



Also, in different cases it is useful to know the result of calculations of the following 
quantity 



Qn{r) 



dx 



r(-r) 



x 



- l - T -TJx) . 



(C.18) 



Taking into account Eq.( |U.16| ), after integration over x and some manipulation we 
obtain for Eq. fldTgp 



r(d/2-r) r(l+n) ^ r(r + A;) 



r(r) r(n + d/2) ^ 2 r(A; + d/2) 

2 2s-k Y(k + n- s + d/2) 



(C.19) 



V(-D* 

r(n-s + l)r 2 (A;-s + l)r(2s-A; + l) 
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Eqs. flU.16l) and ( |U.18| ) allow one to calculate any potential matrix elements. In par- 
ticular, using Eq. flU.l6D , we have very simple expressions for the power-law potential 

2n 

(n r \ : q 2 : \n r ) = — (U.20) 
to 

(n r \ : q 4 : \n r ) = — ^ [d + 6n r - 4] , 

r , i 2n r (n r — 1) r , 
n r : q 6 : \n r ) = V r „ [3d + 10n r - 8 . 

OJ 6 
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pa 


i 


|24 


ORM 


9 = 


0.976562 


e = o 


0.556764 


0.556767 


0.557 


0.5580 






£=l 












1.9375 


9 = 


4.0 


£ = 


2.79573 


2.79575 


2.796 


2.7972 






£=l 












5.65135 


9 = 


62.5 


£ = 


24.79 


24.856 


24. 


86 


24.8569 






£=l 












39.4096 


9 = 


100 


£ = 


34.75 




34.9 


34.9049 






£=l 












54.4364 


9 = 


500 


£ = 


107.0 


108.366 


108 


.4 


108.3854 






£=l 












162.898 


9 = 


1000 


£ = 


174.8 








174.8698 






£=l 












260.295 


9 - 


oo 


£ = 


C = 1.764 


1.85576 


1.856 


1.8559 






£=l 












0.668 






£ = 3 












0.009 



Table 1. The energy of the ground state and orbital excitations 
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c 


£ = Q 


£= 1 


1 = 2 


P 


—E 


P 


—E 


P 


—E 


0.001 


1.0 


0.251 


1.0 


0.06350 


1.0 


0.02877 






(0.251) 




(0.06350) 




(0.02877) 


0.005 


1.0 


0.25497 


1.0 


0.06740 


1.0 


0.032525 






(0.25496) 




(0.06738) 




(0.032520) 


0.01 


1.0 


0.25985 


1.0 


0.07204 


0.99 


0.036813 






(0.25895) 




(0.07202) 




(0.036810) 


0.05 


1.0 


0.29650 


0.99 


0.10219 


0.91 


0.060521 






(0.29650) 




(0.10242 




(0.060520) 


0.1 


0.99 


0.33694 


0.91 


0.12931 


0.88 


0.077046 






(0.33694) 




(0.12931) 




(0.077050) 


0.2 


0.97 


0.404244 


0.88 


0.164292 


0.88 


0.093731 






(0.40424) 




(0.16429) 




(0.093730) 


0.5 


0.91 


0.54246 


0.87 


0.211998 


0.93 


0.107695 






(0.54243) 




(0.21203) 




(0.10770) 


1.0 


0.89 


0.67482 


0.89 


0.236543 


.98 


0.110671 






(0.67482) 




(0.23665) 




(0.11067) 


2.0 


0.89 


0.80554 


0.97 


0.24709 


1.00 


0.111084 






(0.80566) 




(0.24711) 




(0.11108) 


10.0 


0.90 


0.97223 


1.00 


0.249996 


1.00 


0.111111 






(0.97424) 




(0.24998) 




(0.11111) 



Table 2. Energy of ESCP for B = 1 as a function of the parameter c. Result of 
calculations [f| are presented in brackets. 
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£ = 


e= i 


£ = 2 


c 


P 


—E 


p 


—E 


P 


—E 


0.001 


1.0 


0.00168 


1.0 


0.00160 


1.0 


0.001538 






(0.00178) 




(0.00168) 




(0.00160) 


0.005 


0.96 


0.00807 


0.87 


0.00735 


0.88 


0.006746 






(0.00814) 




(0.00739) 




(0.00678) 


0.01 


0.88 


0.015264 


0.72 


0.013575 


0.755 


0.012125 






(0.01540) 




(0.01360) 




(0.01215) 


0.05 


0.69 


0.063234 


0.655 


0.050177 


0.695 


0.040366 






(0.06325 




(0.05019) 




(0.04039) 


0.1 


0.63 


0.11136 


0.67 


0.081862 


0.71 


0.060999 






(0.11130) 




(0.08188) 




(0.06100) 


0.2 


0.64 


0.1885 


0.705 


0.12482 


0.77 


0.083568 






(0.18847) 




(0.12486) 




(0.08360) 


0.5 


0.69 


0.34885 


0.755 


0.18899 


0.88 


0.104978 






(0.34891) 




(0.18915) 




(0.10502) 


1.0 


0.72 


0.51243 


0.86 


0.226511 


.965 


0.110258 






(0.51246) 




(0.22668) 




(0.11027) 


2.0 


0.7863 


0.6877 


0.94 


0.244485 


0.995 


0.111057 






(0.68856) 




(0.24458) 




(0.11106) 


10.0 


0.93 


0.95111 


1.00 


0.24999 


1.00 


0.11111 






(0.95197) 




(0.24997) 




(0.11111) 


Table 3. Similar to Table. 2 for B = 2. 




£ = 


£ 


= 1 


£ 


= 2 


c 


P 


—E 


P 


—E 


P 


—E 


0.005 


0.815 


0.0033 


0.495 


0.0034 


0.455 


0.00345 






(0.00359) 




(0.00355) 




(0.00351) 


0.01 


0.50 


0.0062 


0.465 


0.0065 


0.655 


0.0066 






(0.00700) 




(0.00694) 




(0.00681) 


0.05 


0.50 


0.0311 


0.465 


0.0305 


0.53 


0.02835 






(0.03235) 




(0.03090) 




(0.02840) 


0.1 


0.51 


0.0605 


0.415 


0.0555 


0.535 


0.04780 






(0.06088) 




(0.05570) 




(0.04782) 


0.2 


0.52 


0.110 


0.48 


0.0943 


0.625 


0.07261 






(0.11160) 




(0.09451) 




(0.07260) 


0.5 


0.52 


0.233 


0.64 


0.16370 


0.81 


0.10098 






(0.23359) 




(0.16380) 




(0.10098) 


1.0 


0.53 


0.3790 


0.77 


0.21265 


.94 


0.10951 






(0.37984) 




(0.21279) 




(0.10954) 


2.0 


0.65 


0.56403 


0.90 


0.2402 


0.995 


0.11104 






(0.56407) 




(0.24029) 




(0.11101) 


10.0 


0.94 


0.91317 


1.00 


0.250 


0.667 


0.11111 






(0.91540) 




(0.24993) 




(0.11111) 



Table 4. Similar to Table.2 for B = 4. 
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